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Abstract 

We study how the degree of ordering depends on the strength of the thermal and quantum 
fluctuations in frustrated systems by investigating the correlation function of the order parameter. 
Concretely, we compare the equilibrium spin correlation function in a frustrated lattice which 
exhibits a non-monotonic temperature dependence (reentrant type dependence) with that in the 
ground state as a function of the transverse field that causes the quantum fluctuation. We find the 
correlation function in the ground state also shows a non-monotonic dependence on the strength of 
the transverse field. We also study the real-time dynamics of the spin correlation function under a 
time-dependent field. After sudden decrease of the temperature, we found non-monotonic changes 
of the correlation function reflecting the static temperature dependence, which indicates that an 
effective temperature of the system changes gradually. For the quantum system, we study the 
dependence of changes of the correlation function on the sweeping speed of the transverse field. 
Contrary to the classical case, the correlation function varies little in a rapid change of the field, 
though it shows a non-monotonic change when we sweep the field slowly. 
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I. INTRODUCTION 



Static and dynamic properties of frustrated systems have been studied extensively in the 
past several decades [l-4]. Many model materials have been developed and theoretical stud- 
ies of frustrated systems have been of increasing significance |5|-l8|. In frustrated systems, 
there are highly degenerated ground states and thermal and/or quantum fluctuations gen- 
erate a peculiar density of states. Frustration causes various interesting phenomena such as 
the so-called "order by disorder" and "reentrant phase transition", etc. Fluctuations pre- 
vent the system from ordering in unfrustrated systems, however, in some frustrated systems, 
fluctuations stabilize some ordered structures due to a kind of entropy effect. The ordering 
phenomena due to fluctuations are called "order by disorder" There have been various 



examples of the order by disorder phenomena and also reentrant phase transitions |14N23|. 

Frustration plays an important role in dynamic properties as well as static properties. 
Since there are many degenerate states in frustrated systems, relaxation processes of physical 



quantities often show characteristic features. It is well known that slow re 



random systems such as spin glasses 



2414271] and diluted Ising models (28|-l30|. We also found 



axation appears in 



that slow relaxation processes take place also in non-random systems. We have pointed out 
that frustration causes a stabilization of the spin state by a kind of screening effect, and 
a very slow dynamics appears. We recently proposed a mechanism of slow relaxation in a 
system without an energy barrier for the domain wall 31 1. 

Recently, we have also studied the relation between the temperature dependence of 
static quantities and their dynamics after a sudden change of the temperature in some 
frustrated systems where the temperature dependence of the spin correlation function is 
non-monotonic. In that study, we found that the non-monotonic relaxation of the correla- 
tion function indicates a picture of relaxation of an effective temperature of the system 32 |. 

Quantum fluctuations also cause some ordering structure in frustrated systems as well as 
thermal fluctuations 33|]. The so-called quantum dimer model is the simplest model that 
exhibits the order by quantum disorder. This system has a rich phase diagram as a function 
of the on-site potential and the kinetic energy. When the on-site potential and the kinetic 
energy are equal, the model corresponds to a Rokhsar-Kivelson point 3J], while the model 
corresponds to an Ising antiferromagnet with the transverse field on its dual lattice, when 
the on-site potential is zero. Recently, "order by quantum disorder" was also studied in a 
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transverse Ising model on a fully frustrated lattice from a viewpoint of adiabatic quantum 



annealing 



35|. 



The nature of the dynamics in quantum systems with time-dependent external fields has 



also been studied extensively 
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3H40|. It is an important issue to investigate the real-time 
dynamics in quantum system for control of the quantum state by external fields. Recently, 
quantum annealing or quantum adiabatic evolution have been studied from a viewpoint of 



quantum dynamics 



4l|-|44l|. They are the methods to obtain the ground state in complex 



systems. Now, the quantum annealing method is adopted on a wide scale, for example, 



clustering problem and Bayes in 
and information engineering 45 



erence which are important topics in information science 



46|. 



The purpose of the present study is to understand the similarity and difference between 
the effects of the thermal fluctuations and the quantum fluctuations in frustrated systems. 
In this paper, we study the dynamic properties of a frustrated Ising spin system with a 
time-dependent transverse field. In Section 2, we introduce a frustrated Ising system with 
decorated bonds, and we review equilibrium and dynamical properties of this classical sys- 
tem. In Section 3, we study the ground-state properties of this model with a transverse 
field, and we investigate the dynamical nature of this model with a time-dependent trans- 
verse field. We also consider the microscopic mechanism of the non-monotonic dynamics of 
the spin correlation function. In Section 4, we summarize the present study. 



II. CLASSICAL DECORATED BOND SYSTEM 



A. Equilibrium Properties 

First we briefly review static properties of the model depicted in Fig. [H The circles and 
the triangles denote the system spins o"i and 02 and the decoration spins o"j (z = 3, ■ ■ ■ + 
respectively, where is the number of the decoration spins. Both of the system spins and 
the decoration spins are 5 = 1/2 Ising spins. The Hamiltonian of this system is given by 

= ^(^1(^2 - J ^ (o-i + (Ja) o-j, (1) 

where Jo and J are positive. The direct bond Jo tends to have an antiferromagnetic corre- 
lation while the interactions through the decoration spin J tend to have ferromagnetic one. 
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We will study effects of the competition between them. If we set Jq = N^J, the interactions 
cancel out each other at zero temperature. Because the interaction through the decoration 
spin is weakened at finite temperatures as we will see below, the spin correlation is antifer- 
romagnetic at finite temperatures in this case. In this paper we focus the case where the 
spin correlation function changes the sign at a temperature, and therefore we study the case 
Jo < NdJ. Concretely, we adopt the case Jo = NdJ/2. The consequences obtained in the 
present paper hold qualitatively for any choice of the ratio Jq/N^J as long as it is less than 
unity. We take J as the energy unit from now on. Figure [2] shows the equilibrium spin 
correlation function between the system spins ai and (J2 (hereafter, we call it the "system 
correlation function" and denote it by C) as a function of the temperature in the cases of 
Nd = 1, 2, 4, 6, and 8. In Fig. [21 Tj^in and To represent the temperatures where C takes the 
minimum value and C = 0, respectively. In the present case, Tmin ~ 3.64 and To ~ 1.64. 

The system correlation function C behaves non-monotonically as a function of tempera- 
ture. This non-monotonic behavior comes from a kind of entropy effect. At T = 0, all the 
spins align in the same direction, because energetically it is the most favorable state. At 
finite temperatures the decoration spins can flip due to the thermal fluctuations. When each 
decoration spin has the values ±1 randomly, the ferromagnetic interactions depicted by the 
solid lines in Fig. [1] are weakened, and the direct antiferromagnetic interaction becomes dom- 
inant. The states, in which the decoration spins align randomly, are entropically favorable. 
Therefore, at a high temperature, C becomes negative due to the entropy effect. This is why 



0-1 




0-2 



FIG. 1: Decorated bond system. The solid and the dotted lines represent ferromagnetic inter- 
actions —J and an antiferromagnetic interaction Jo, respectively. The circles and the triangles 
denote the system spins, cJi and cj2, and the decoration spins cjj (i = 3, • • • , A'd + 2), respectively. 
This figure shows the case of A^a = 4. 
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C behaves non-monotonically as a function of the temperature. Following Refs. 14|-|23| , the 
effective Hamiltonian Hcfr = — Jcfro'io'2 is defined by 

^ ... e-'^^ = (2cosh2/3J)^e-^^^«. (2) 

o-3=±l (TjVj^2=±l 

The reduced coupling of the system spins is obtained as 

i^eff = /3 Jeff = ^ log cosh (2/3 J) - J. (3) 

The temperature dependence of C can be calculated analytically by tracing out the degree 
of freedom of the decoration spins: 

C = (0-10-2) = tanhi^'eff- (4) 

Here it should be noted that the effective interaction Eq. (jl]) is proportional to N^, and 
thus the temperature at which C = is the same for all A^d and that temperature of the 
minimum C as well. 



B. Kinetic dynamics 



In the previous section, we showed the equilibrium properties of the decorated bond 
system. In this subsection, we study relaxation processes of C after the temperature is 



To T, 





N,= l 




Nd = 2 : 




Nd = 4 




Nd = 6 _ 




Nd = 8 — 
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FIG. 2: (Color online) The system correlation function C given by Eq. (j4]) between the system 
spins cJi and (T2 as a function of the temperature. In this case, Tmin — 3.64 and Tq ~ 1.64. 
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decreased suddenly. We adopt the Glauber Ising model for the time evolution: 

dt 

i 
i 

where P (ai, ■ ■ ■ , Uj, ■ ■ ■ , t) = P ({ci} ; t) and Wo-._j._o-. represent the probability distri- 

bution at time t and the transition probability from cr^ to — Uj, respectively. The transition 
probability Wai^-ai is given by 



_ -t eg yul, - ■ ■ , —Uj, • • • , /„x 



Peq ((Ti, ■ ■ ■ , (Tj, ■ ■ ■ , (T7V^+2) + -Pcq ((Tl , ■ ■ ■ , — (Tj, ■ ■ ■ , CrAr^+2) 

where Poq (ai, ■ ■ ■ , (Xj, ■ ■ ■ , = -Pcq ({cj}) denotes the equilibrium probability distri- 

bution at the temperature T. 



In Ref . |32[ , we studied the dynamics of C after a sudden change of the temperature from 
a finite temperature Ti to a finite temperature T2. In this paper, we change the temperature 
suddenly from Ti = 00 to a finite temperature T2. The dynamics of C is given by 

C{t) = J2PiW^ht)a^a2, (7) 

where P {{(Ji} ; t) is the probability of the configuration {ai} at the time t. 

From the equilibrium properties, it is trivial that C behaves non-monotonically, when 
we decrease the temperature slowly enough. However, it is not clear whether C behaves 
monotonically or not when the temperature is changed suddenly. 

The initial condition is set in the equilibrium state at Ti = 00, in other words, the initial 
condition is the uniform distribution: P({(Tj};t = 0) = 1/2^'^"'"^. We suddenly decrease 
the temperature to a finite temperature T2. Figure [3] shows the time evolutions of C for 
iVd = 2,4,6,8, and 10 in the cases of T2 = 7.29, 3.64(~ Trnm), 2.35, 1.64(~ Tq), 1, and 0.5. 
As we saw in Fig. |2l at T2 = 7.29, 3.64, and 2.35, the equilibrium value of C is negative. At 
T2 = 1.64, the equilibrium value of C is zero, and at T2 = 1.0 and 0.5, it is positive. 

In the case of T2 = 7.29 > Tmin, C monotonically decreases toward the equilibrium value. 
On the other hand, non-monotonic relaxations appear in the cases of T2 = 3.64, 2.35, 1.64, 1, 
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(c)T2 = 2.35 
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FIG. 3: (Color online) Dynamics of C for T2 = 7.29, 3.64(~ Tmin), 2.35, 1.64(~ Tq), 1, and 0.5. 
The initial condition is set to be the equilibrium probability distribution at T = 00. The lines have 
the same interpretations as in Fig. [2j Note the different scales on the axes in the different parts. 



and 0.5 which are lower than Tmin. Such a nature was pointed out in Ref. 32|. Although the 
temperature T is suddenly changed to a low temperature T2, the property of the system, i.e., 
the correlation function, shows a gradual change which seems to follow the gradual decrease 
of the temperature. Therefore, we regard that a kind of "effective temperature" of the 
system decreases gradually. This non-monotonic behavior takes place as long as T2 < Tmin. 
This fact indicates that such non-monotonic relaxation does not depend on the equilibrium 
value of C at the final temperature T = T2. It should be noted that as increases, the 
time evolution shows stronger non-monotonicity. 
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III. DECORATED BOND SYSTEM WITH TRANSVERSE FIELD 



In the previous section, we introduced the decorated Ising model and studied its equihb- 
rium and dynamical properties. In this section, we study the ground-state properties of this 
model with a transverse field, and the dynamical nature of this model with a time-dependent 
transverse field. 



= -^Jalal -JY,K + ol) (9) 



A. Ground-State Properties 

The Hamiltonian is given by 

n{a)=an^\{\-a)U^ (0 < a < 1), (8) 

where 

= - 

2 

i=3 

and 

i=l 

Here, "He and "Hq represent the classical part and the quantum part of the total Hamiltonian, 
respectively. The operators a? and af are the Pauli matrices. 

Hereafter, we call a the "quantum parameter". The total Hamiltonian (Eq. (|8])) with a = 1 
corresponds to the completely classical Hamiltonian. The total Hamiltonian possesses parity 
symmetry. Now we make use of the all-spin-flip operator: 

V=l[a: (12) 

i=l 

to study the parity symmetry. Total Hamiltonian and the all-spin-fiip operator commute: 

[n{a),V] = 0, (forVa) (13) 
Symmetric wave function and antisymmetric wave function are defined as follows: 

l*s) = E'w«.(l^)+^k))> (14) 

l<^'as)=Ew«.(l^)-^k)), (15) 
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FIG. 4: (Color online) Eigenenergies as functions of a for = 1, 2, and 4. Solid and dotted 
curves denote the eigenenergy for symmetric wavefunctions and for antisymmetric wavefunctions, 
respectively. 



spin configurations fixing cti = +1. |$s) and |$as) represent symmetric and antisymmetric 
wave functions, respectively. Because V\^s) = I'^'s) and "Pj^as) = — I'^'as) are satisfied, 
the total Hamiltonian can be block-diagonalized according to the symmetry. Note that the 
ground state of "H (a) with arbitrary a is a symmetric wave function. Figure H] shows the 
eigenenergies as functions of the quantum parameter a for N^^ = 1,2, and 4. The sohd and 
the dotted curves in Fig. H] denote the eigenenergy values for symmetric and antisymmetric 
wavefunctions, respectively. 

We consider the system correlation function in the ground state: Cgs (a) = 
(\l/gs (a)| crfcrf |\l/gs (a)), where |\I^gs(a)) denotes the ground state of the Hamiltonian 
(Eq. (jH])) with the quantum parameter a. Figure [5] shows the system correlation func- 
tion in the ground state as a function of the quantum parameter a. The system correlation 
function of the ground state behaves non-monotonically as well as the equilibrium value 
at the finite temperature as shown in Fig. [2j This fact indicates a similarity between the 
thermal fluctuation and the quantum fluctuation, although some details are different, e.g. 
the points where the system correlation function Cgs becomes zero and the minimum values 
depend on the number of the decoration spins in the quantum case, whereas they do not 
depend on the number of the decoration spins in the classical case. 

To consider the microscopic mechanism of this non-monotonic behavior, we calculate the 
probability of the classical bases 



where \a) denote a state of a spin configuration, and Yl'{a} denotes summation over all the 




(16) 
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where \i) = |+, +,+,■■ ■ , +) and |— ,+,+,■■ ■ , +) etc. Figure [6] shows the probabihty dis- 
tribution as a function of the quantum parameter. The sohd and the dotted hues in Fig. [6] 
denote P{i;a) for the states in which o"fcr| = +1, and afcr| = —1 states, respectively. In 
the case of N^i = 1, though there are eight classical configurations: (ai, (72, a^) = (+, +, +), 
(-,+,+), (+,-,+), (-,-,+), (+,+,-), (-,+,-), (+,-,-), and (-,-,-), only three 
configurations (+, +, +), (+, +, — ), and (+, — , — ) give different values because of the sym- 
metry. Thus we see only three lines in the left panel of Fig. |6l By the same reason, we see 
only five lines and eight lines in the middle and the right panels of Fig. [61 



B. Real-Time Dynamics 

In the previous section, we studied the ground-state properties of the decorated bond 
system with the transverse field. The system correlation function behaves non-monotonically 
as a function of the quantum parameter a. Now, we consider the real-time dynamics of 
the system correlation function in the quantum case by the time-dependent Schrodinger 
equation. The system correlation function C{t) is defined by 



C{t) = {^{t)\ala^,\^{t)) 



(17) 



where \ip (t)) denotes the wavefunction at time t. Now we consider the time-dependent 
Hamiltonian expressed by 



n (t) = -Uc 

T 



n 



(18) 
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FIG. 5: (Color online) The system correlation function of the ground state Cgs in the ground state 
as a function of the quantum parameter a. 
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FIG. 6: (Color online) Adiabatic limit of the probability distribution function P {i; a) for iVj = 1, 2, 
and 4 (from left to right). The solid and the dotted lines denote the probability of the states where 
the system correlation functions are positive and negative, respectively. Some of P{i;a)s give the 
same curves and they are overlapped each other. 

where represents the sweeping speed. Here, t/r corresponds to a in the previous section. 
The initial condition is set to be the ground state of "H (t = 0) = "Hq such as 

|^(t = 0)) = K,---,^), (19) 

where |— >) = (|t) + \i))/V^- Because this ground state is a symmetric wavefunction, the 
wavefunction is also symmetric after the sweeping. If the sweeping speed is slow enough 
[i.e. the adiabatic limit), the system correlation function behaves as given in Fig. [61 

Figure [7] (a)-(e) shows the real-time dynamics of the system correlation function C{t) 
in the cases of = 1, 2, 4, 6, and 8 for several values of r. For large values of r, C{t) 
is non-monotonic as a function of time, while it changes little when r is small. The final 
value C{t) is non-monotonic as a function of r. For the largest value of r, C(t = r) ^ 1 
which is the value in the ground, while it is nearly zero for the small value r = 0.1. In 
between, C(r) takes negative values. The change becomes significant when increases. 
We summarize these dependence in Fig. [7](f). Here we see that, for small values of r, C(r) 
moves to a negative value monotonically as r increases, and that C (r) has a minimum point 
at an intermediate value of r. 

The negative C(r) is considered to be attributed to the negative value of Cgs at an inter- 
mediate value of a. And thus, we again find a gradual change of the effective temperature 
of the system as we saw in the case of temperature changing protocol of the classical system 
as stated in Section fll B[ However, if we change a fast (small r), we find monotonic changes 
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FIG. 7: (Color online) Figures (a)-(e) give the real time dynamics of the system correlation 
function for = 1, 2, 4, 6, and 8, respectively, (f) The final values of the system correlation 
function as a function of r (semilog scale) . 



which are very different from the classical case. This indicates a difference between thermal 
fluctuation and quantum fluctuation effect. 

Next we consider the microscopic mechanism of this non-monotonic dynamics. We 
calculate the real-time dynamics of the probability P{i) of the classical bases [i = 
(+,+,+),■■■ ,(—,—,—)) for = 1,2, and 4 for r = 1 and 10. The results are shown 
in Fig. IHl There are curves, but many of them are the same and overlap. As r increases, 
P{i)s approach the adiabatic ones shown in Fig. [61 

If we consider probability P_ to have crj^o"! = —1 which is the sum of 
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FIG. 8: (Color online) Probability P{i) for classical basis states i in case of iVd = 1,2 and 4 for 
r = 1 and 10. The curves have the same interpretations as in Fig. [6l 



P(— , +, +),P(— , — , +),P(+, +, — ),and P(+,— ,— ), it is zero in the adiabatic limit. How- 
ever, in the short time, for r = 1, P_ remains nonzero. For r = 10, P_ is close to zero. 

We also calculate the dynamic behavior of the probability of the eigenstate {/} at several 
times t/r. Here, the eigenstates are labeled in the order of eigenenergy (/ = 1,2,- ■■ ,8). 
Because the wavefunction after changing the Hamiltonian is symmetric, it is only necessary 
to consider the eigenenergy of symmetric wave functions. Figure [9] shows the dynamical 
behavior of the probability at the state {/} in cases of A^d = 1, 2, and 4 for r = 1 and 10. 
For T = 1, the probability of the first excited level is larger than the probability of the ground 
state. On the other hand, for r = 10, the probability distribution decreases monotonically as 
a function of the eigenenergy. As r increases, the probability of the ground state approaches 
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FIG. 9: (Color online) The probability P(Z) of eigenstates (/ 
indicates the ground state. 
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IV. CONCLUSION 



We studied a correlation function in a frustrated quantum spin system with a transverse 
field. In the ground state, the system correlation function is non-monotonic as a function of 
the quantum parameter a, as is the equilibrium value of the system correlation function as a 
function of temperature in the classical case. We also considered the dynamics of the system 
correlation function of the system with a time-dependent quantum parameter. When we 
increased the quantum parameter a fast, the wave function changed little from the initial 
wave function, and the dynamics is monotonic. As the sweeping speed of the quantum 
parameter a decreases, the dynamics of the system correlation function approaches the 
adiabatic limit, and the non-monotonic relaxation appears. In the classical system, however. 
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we found non-monotonic relaxation of the system correlation function when we decreased 
the temperature suddenly. This fact shows the difference between the effects of the thermal 
fluctuation and the quantum fluctuations. In frustrated systems, correlation function often 
behaves non-monotonic due to a peculiar density of states. We expect that the mechanism 
which is discussed in this manuscript appears in real frustrated systems. 
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